Emergence of half-metallicity in suspended NiO chains 
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Contrary to the antiferromagnetic and insulating character of bulk NiO, one-dimensional chains 
of this material can become half-metallic due to the lower coordination of their atoms. Here we 
present ab initio electronic structure and quantum transport calculations of ideal infinitely long 
NiO chains and of more realistic short ones suspended between Ni electrodes. While infinite chains 
are insulating, short suspended chains are half-metallic minority-spin conductors which display very 
large magnetoresistance and a spin-valve behaviour controlled by a single atom. 
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In going from bulk to lower dimensions material prop- 
erties often change drastically. A recent example is that 
of interfaces between different insulators which can ac- 
tually be metallicQ. Even more recently, it has been 
predicted theoretically that certain oxygen surfaces of 
some insulating ceramic oxides can exhibit magnetism 
and half-metallicity Q. The ultimate limit in this re- 
spect can be found in atomic chains formed in metallic 
nanocontacts which allow to study the transport prop- 
erties of one-dimensional systems of atomic size (3}- Due 
to the lower coordination of the atoms the properties of 
metallic atomic chains formed in nanocontacts can be re- 
markably different from those in the bulk. For example, 
Pt nanocontacts can exhibit magnetic order when atomic 
chains are formed 0- More complex one-dimensional sys- 
tems like carbon-cobalt atomic chains^ or organometal- 
lic benzene-vanadium wiresQ have even been predicted 
to be half-metallic conductors. 

However, not all metals form atomic chains in nanocon- 
tacts, although recently, it has been found that the pres- 
ence of oxygen favours their formation]?]. For exam- 
ple, experiments with Ni nanocontacts ja 111 have never 
shown evidence of chain formation. Nevertheless, the 
presence of oxygen in the contact region could possibly 
lead to the formation of NiO chains. In this context it 
has also been proposed that the rather moderate mag- 
netoresistive properties of pure Ni nanocontacts [ffj, |ll| 
could be enhanced considerably by the presence of oxy- 
gen adsorbates on the surface of the Ni electrodes|l2fl. 
On the other hand, bulk NiO is a common example of 
a correlated insulator with antiferromagnetic (AF) order 
(see, e.g., Ref. El El, 

which remains insulating even 
above the Neel temperature when the AF order is lost. 
Thus it is not at all obvious whether or not oxidized Ni 
nanocontacts or NiO chains should be conductors. 

In this Letter we investigate the electronic and mag- 
netic structure and the transport properties of one- 
dimensional NiO chains, both idealized infinite ones and 
more realistic short ones suspended between Ni nanocon- 
tacts. Anticipating our most important results our ab 
initio quantum transport calculations show that short 



NiO chains suspended between Ni nanocontacts can be- 
come half-metallic conductors, i.e., carry an almost 100% 
spin-polarized current. This result holds true even for a 
single O atom in between Ni electrodes. Consequently, 
for antiparallel alignment of the electrode magnetizations 
the transport through the contact is strongly suppressed 
resulting in very large MR [difference in resistance be- 
tween antiparallel and parallel alignment of the mag- 
netizations of the electrodes normalized either to the 
higher resistance value (MRi) or to the lower one (MR2)]: 
MRi « 90% and MR 2 » 700%, respectively. 

The electronic structure of both idealized infinitely 
long NiO chains and the more realistic short suspended 
chains have been calculated in the density functional 
theory (DFT) approximation using the hybrid func- 
tional B3LYP|lfij. which combines local and non-local 
exchange. This functional has been successful in de- 
scribing the electronic and magnetic structure of some 
strongly correlated materials like bulk NiO. In particular, 
it describes very well the charge-transfer character, the 
magnitude of the gap and the magnetic moment of NiO, 
and predicts the correct AF order [li). It is important 
to stress here that due to the insufficient cancellation of 
the self-interaction in the local exchange functional of the 
local density approximation (LDA) the occupied narrow 
3d bands are raised in energy. As a result LDA severerly 
underestimates the gap of bulk NiO 01 • The generalized 
gradient approximation (GGA) to the exchange func- 
tional improves somewhat the description of NiO but the 
energy gap of ^2eV is still too small and also the charge- 
transfer character is not captured correctly E3 ■ B3LYP 
corrects the self-interaction error inherent in LDA and 
GGA by mixing the exact non-local Hartree-Fock (HF) 
exchange with the GGA exchange functional. The DFT- 
based alternatives that correct the self-interaction error 
like LDA+U, the self-interaction-corrected LDA (SIC- 
LDA), and the GW approximation lead to similiar results 
as those obtained from B3LYP 14J . 

The electronic structure of infinite NiO chains has been 
calculated with the CRYSTAL03 program package E3 
while for the electronic structure and quantum transport 
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FIG. 1: Schematic one-electron energies of a one-dimensional 
NiO chain in z-direction for minority-spin. To the left and 
right the orbital energies of an individual Ni + cation and O - 
anion in the crystal field of a one-dimensional Ni + 0~ chain 
are shown. In the center the formation of valence and con- 
duction bands by hybridization of Ni 3d and O 2p orbitals is 
shown. 



calculations of the suspended chains we have used our 
ab initio quantum transport package AL AC ANT [Til. Hfljl . 
The infinite chain calculations were done with elaborate 
all-electron basis sets for Ni and O I20L I21I , similiar to 



those employed for reported HF and B3LYP calculations 
of bulk NiO|l4|, but extended with a diffusive sp- function 
in the case of Ni and a d-polarization function in the case 
of O. Using these basis sets we reproduce the B3LYP re- 
sults for bulk NiO We have also employed minimal 
basis sets with effective core pseudo-potentials described 
in earlier work 0] for both Ni and O, and have found 
that the results change very little. Thus we have em- 
ployed these minimal basis sets for the more demanding 
transport calculations. 

The electronic properties of NiO are, to a large ex- 
tend, determined by the atomic scale properties, like the 
crystal field splitting of the Ni and O energy levels and 
the amount of electron charge transfered from Ni to O. 
The latter, in turn, is determined by the interplay be- 
tween Madelung binding energy, the ionization potential 
of Ni, and the electron affinity of O. Due to the lower co- 
ordination and the corresponding decrease in Madelung 
binding energy the electron transfer from Ni to O is less 
favourable in an atomic chain than in bulk (where the 
electron transfer is almost complete resulting in an ionic 
configuration of Ni 2+ 2 ~). The proper starting point 
to discuss the formation of energy bands in the one- 
dimensional NiO chain are therefore the univalent ions 
Ni+ and O". 

In order to understand how the low coordination af- 
fects the atomic properties of the constituting ions we 
have first performed B3LYP calculations of both a single 
Ni + ion and a single O - ion each in the field of point 
charges that mimic the crystal field of a one-dimensional 
chain of univalent Ni and O ions. We find that for both 
ions the spin-doublet state (S = 1/2) minimizes the en- 
ergy as is the case for the free ions. In Fig. \I\we show 
schematically the energy levels of Ni + and O - in the pres- 
ence of the point charges for minority spin only. Inter- 
estingly, for minority spins, the occupied Ni 3di orbitals 
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FIG. 2: (a) Energy per unit cell of infinite NiO chain in 
dependence of lattice spacing a. Black triangles indicate the 
half- metallic state (HM), grey triangles the molecule-like in- 
sulating state with FM order (ML), and black boxes the in- 
sulating state with AF order (AF). (b) Band structure of 
HM state for a lattice spacing of 3.6A. Solid lines indicate 
majority-spin bands and dashed lines indicate minority-spin 
bands, (c) Same as (b) but for ML FM state. 



(d xz and d yz ) of the Ni + ion fall energetically in between 
the occupied and unoccupied O 2pi orbital (j> x and p y ) 
of the O - ion. Thus the 3di and 2pi orbitals can form 
two filled degenerate bonding bands and two degener- 
ate partially filled antibonding bands as indicated in the 
middle part of Fig. ^ The 3c?2 (d xy and d x i_ y i) doublet 
is somewhat above in energy to the 3d\ but cannot hy- 
bridize with the oxygen 2p orbitals. The Ni 3do {d^ z i_ r i) 
and 4s are empty while the O 2po (p z ) and 2s orbitals are 
filled and much lower in energy so that no hybridization 
takes place though symmetry would allow for it. 

On the other hand, for majority-spin electrons (not 
shown) all five Ni 3c? orbitals are filled while the 4s is 
also empty, i.e., the Ni + valence configuration is 3d 9 and 
not 4s 1 3d 8 as for the free Ni + ion. Moreover, all of the 
O 2p orbitals are filled so that the Ni + and 0~ ions can 
only form either completely filled or completely empty 
bands for the majority spin. Thus the ionic picture sug- 
gests that a one-dimensional NiO chain should become a 
half- metallic conductor where only one kind of spin levels 
forms conducting bands. In this case only the minority- 
spin 3di bands would conduct so that, according to the 
proposed classification in Ref. |U it would be of Type 
Ib. 
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Not surprisingly, our calculations for inifinite onc- 
dimensional NiO chains (Fig. 0) show that a univalent 
ionic configuration as an initial guess results in a half- 
metallic state for large separation of the individual chain 
atoms (i.e., large lattice spacing of 5 A) as suggested by 
the ionic picture. However, as can be seen from Fig|^a), 
the half-metallic state is only a metastable state for most 
values of the lattice spacing. This half-metallic state is 
"shadowed" by a second state with FM order and insu- 
lating character. By successively decreasing the lattice 
spacing a of the chain and using the (half- metallic) state 
of the previous step for the initial guess, the half-metallic 
state can be generated also for smaller inter-atomic dis- 
tances which points towards its metastability. Around 
the equilibrium lattice spacing (a ~ 3.4A) the ground 
state of the chain has AF order and is of insulating char- 
acter with a substantial gap of ~ 4eV like in bulk. When 
stretched out of equilibrium the FM state and the AF 
state become comparable in energy until finally, at a lat- 
tice spacing of ~ 4.2A, the FM state becomes the ground 
state. 

The band structure diagram in Fig. |2fc) shows that 
the metastable state with FM order corresponds indeed 
to the half-metallic state suggested by the ionic pic- 
ture: The half-filled doubly-degenerate conduction band 
is formed by minority-spin Ni 3d% orbitals hybridized 
with O 2pi orbitals while the Ni 3c?2 orbitals do not hy- 
bridize with O 2p orbitals and thus form a flat valence 
band. The lowest-lying empty band is formed by the 
minority-spin Ni 3c?o orbital which is slightly hybridized 
with the Ni 4s orbital. On the other hand the stable state 
with FM order and insulating character (see band struc- 
ture in Fig. Htb)) actually corresponds to the ground 
state of the NiO molecule which is a 3 £~ 1 state[23j. 
The main difference with the half-metallic state is that 
now the non-degenerate Ni 3drj and 4s orbitals form a 
minority-spin valence band while the minority-spin dou- 
bly degenerate half-filled antibonding band composed of 
Ni 3di and O 2pi bands is emptied and a substantial gap 
of ~ 3eV opens. Thus the infinite chain behaves like an 
insulator for reasonable values of the chain stretching. 

Atomic chains formed in break junctions have a finite 
length and are suspended between electrodes. It is well 
known that the contact between the atomic chain and 
the electrode tip will have considerable effect on the elec- 
tronic structure of the chain, especially when <i-orbitals 
are involved like is the case here|lfj|. We have thus calcu- 
lated the electronic structure and transport properties of 
both a single oxygen atom and a O-Ni-0 chain bridging 
the two tips of a Ni nanocontact as shown in the insets of 
the right panels of Fig. |3]and FigQ] In the case of the sin- 
gle oxygen atom (Fig. |3J) the electron transport is almost 
100% spin-polarized around the Fermi level for parallel 
(P) alignment of the magnetizations of the two Ni elec- 
trodes. Moreover, an orbital eigenchannel analysis jiij 
reveals that the transport is mainly due to two almost 
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FIG. 3: Transmission per spin channel in the case of P (left) 
and AP (right) alignment of the electrode magnetizations for 
NiO chain consisting of one oxygen atom bridging the two Ni 
tip atoms of the Ni electrodes as shown on the inset in the 
right panel. The separation of the two Ni tip atoms is 3.6A. 



perfectly transmitting minority-spin channels composed 
of Ni 3di and O 2pi orbitals, i.e. they correspond to 
the doubly-degenerate conduction band of the metastable 
half-metallic state in the perfect chain. Thus the half- 
metallic state which was suppressed in the idealized case 
of the infinite chain emerges in the more realistic situa- 
tion of a short suspended chain. We can understand this 
phenomenon in terms of the orbital blocking mechanism 
proposed earlier in the context of Ni nanocontacts [1(1, l24j . 
The highest minority-spin valence band of the insulating 
state with FM order in the infinitely long chain has a 
strong contribtution from the Ni 3rio orbital which is not 
"compatible" with the geometry of the pyramid shaped 
Ni contacts, so that this band is blocked and thus cannot 
be occupied. Instead, the doubly-degenerate band com- 
posed of Ni 3di orbitals hybridized with O 2pi orbitals is 
partially filled resulting in the half-metallic state which 
in the perfect chain is only metastable. Thus the orbital 
blocking by the contacts actually turns the chain into a 
half-metallic conductor. Consequently, the conductance 
is strongly suppressed in the case of antiparallel (AP) 
alignment of the magnetizations of the Ni electrodes as 
can be seen from the right panel of Fig. [3] and the MR 
becomes very large: MRi « 90% and MR 2 » 700%. 

Varying the distance d between the Ni tip atoms leads 
to similiar results as those shown in Fig. [3J For the P 
case the current through the chain is almost 100% spin- 
polarized with two open minority channels composed of 
Ni 3g?i and O 2pi orbitals, while for the AP case it is 
strongly suppressed, resulting in very high MR values 
between 80% and 90% for MRi for d between 3.0A and 
5.0A. Geometry relaxations for different values of the tip- 
tip distance show that for small distances the oxygen 
atom goes into a zigzag position. The bonding angle 
decreases with increasing distance until it becomes zero 
at 3.6A. Finally, the chain breaks for d > 4.8A. Thus the 
scattering is strong for small distances d < 3.6A when 
the bonding angle is appreciable and for large stretching, 
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FIG. 4: Transmission per spin channel for suspended chain 
consisting of a O-Ni-O bridge and the Ni tip atoms shown in 
inset of right panel in the case of FM order (left) and of AF 
order (right). For AF order the magnetization of the Ni atom 
in the center of the chain is reversed with respect to the two 
Ni tip atoms. Distance between a Ni tip atom and the center 
atom is 3.6A. 



sion of the current for AP alignment of the electrodes 
resulting in very large MR values of MRi f» 90% and 
MR 2 w 700%, respectively. This could perhaps explain 
to some extend the very large MR values in Ni nanocon- 
tacts obtained in some experiments^ where oxygen is 
likely to be present. Finally, the O-Ni-0 bridge sus- 
pended between Ni electrodes operates as a single atom 
spin valve where the currentflow is controlled by the mag- 
netization of a single atom. 
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cial support from Grant No.FIS200402356 (MCyT), No. 
GV05-152, and Ramon y Cajal Program (MCyT). J.J.P. 
acknowledges financial support from MAT2005-07369. 



d > 4.2A, resulting in a considerable reduction in the 
conduction of the two minority channels. 

In longer suspended chains the insulating state with 
FM order starts to emerge inside the chain and a away 
from the contacts. As a result the conductance is re- 
duced as can be seen already in the case of the O-Ni-0 
bridge (Fig. 0J). The minority-spin conduction is reduced 
considerably (~ 30%) compared to the case of the single 
oxygen bridge. On the other hand the conduction of the 
majority-spin channel becomes practically zero (< 0.2%). 
This can be understood by the fact that the residual 
majority-spin channel conductance of the single oxygen 
bridge is due to direct hopping of Ni s electrons between 
the electrodes and therefore vanishes when the distance 
between the electrodes is large. The finite minority-spin 
conductance opens up the possibility to an interesting 
phenomenon: When the middle Ni atom reverses its spin, 
the conductance drops to nearly zero (see right panel of 
Fig. 0} since the AF chain is insulating. In other words, 
this system behaves as a single atom spin valve which 
presents an extremely large MR even higher than that 
reported above for the single oxygen bridge: MRi w 99% 
and MR2 « 10, 000%. Apart from controlling the magne- 
tization direction of the central atom by a magnetic field, 
Fig. Ota) suggests that a mechanical control of the spin 
valve (by stretching the chain) would also be possible. 

In conclusion, we have shown that one-dimensional 
infinite NiO chains are insulating and have AF order 
around the equilibrium lattice spacing but can actually 
become a FM insulator when stretched out of equilib- 
rium, in contrast to bulk NiO. In the more realistic case 
of a short NiO chain suspended between Ni nanocontacts 
the chain becomes an almost 100% spin-polarized con- 
ductor. The emergence of almost perfect half-metallicity 
(i.e. almost 100% spin-poalrization of the conduction 
channel) in suspended chains leads to a strong suppres- 
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